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Abstract 

In agreement with the Harris criterion, asymptotic critical exponents of three- 
dimensional (3d) Heisenberg-like magnets are not influenced by weak quenched 
dilution of non-magnetic component. However, often in the experimental studies 
of corresponding systems concentration- and temperature-dependent exponents are 
found with values differing from those of the 3d Heisenberg model. 

In our study, we use the field-theoretical renormalization group approach to ex- 
plain this observation and to calculate the effective critical exponents of weakly 
diluted quenched Heisenberg-like magnet. Being non-universal, these exponents 
change with distance to the critical point T c as observed experimentally. In the 
asymptotic limit (at T c ) they equal to the critical exponents of the pure 3d Heisen- 
berg magnet as predicted by the Harris criterion. 
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1 Introduction 



Relevance of structural disorder for the critical behaviour remains to be an 
important problem of modern condensed matter physics. Even a weak disor- 
der may change drastically the behaviour near the critical point and in this 
respect may be related to the global characteristics of a physical system, such 
as the space dimension, order parameter symmetry and the origin of interpar- 
ticle interaction. In this paper, we are going to discuss some peculiarities of a 
paramagnetic-ferromagnetic phase transition in magnets, where the random- 
ness of structure has the form of substitutional random-site or random-bond 
quenched disorder. Solid solutions of magnets with small concentration of non- 
magnetic component as well as amorphous magnets with large relaxation times 
may serve as an example of such systems. 

Intuitively, it is clear that for a weak enough disorder the ferromagnetic phase 
persists in such systems. Obviously, intuition fails to predict whether the crit- 
ical exponents characterizing phase transition into ferromagnetic state will 
differ in a disordered system and in a "pure" one. The answer here is given 
by the Harris criterion [1] which states that the critical exponents of the dis- 
ordered system are changed only if the heat capacity critical exponent of a 
pure system is positive, otherwise the critical exponents of a disordered system 
coincide with those of a pure one. Returning to d — 3 dimensional magnets 
with 0{m) symmetric spontaneous magnetization one is lead to the conclu- 
sion, that here only the critical exponents of uniaxial magnets described by 
the d = 3 Ising model (m = 1) are the subject of influence by weak quenched 
disorder. Indeed, the heat capacity diverges a = 0.109 ± 0.004 > [2] for 
m — 1, whereas it does not diverge for the easy-plane and Heisenberg-like 
magnets: a = -0.011 ± 0.004 and a = -0.122 ± 0.010 for m = 2 and m = 3, 
respectively [2]. 

Note however that the Harris criterion tells about the scaling behaviour at the 
critical point T c . In other words it predicts (possible) changes in the asymptotic 
values of the critical exponents defined at T c . In real situations one often 
deals with the effective critical exponents governing scaling when T c still is 
not reached [3]. These are non- universal. As far as in our study of particular 
interest will be the isothermal magnetic susceptibility xt let us define the 
corresponding effective exponent by [3]: 



In the limit T — > T c the effective exponent coincides with the asymptotic one 

7eff = 7- 



7cff(r) 



dlnx(r) 
dlnr 



with r=\T-T c \/T c . 
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Already in the first experimental studies of weakly diluted uniaxial (Ising-like) 
d = 3 random magnets [4] the asymptotic values of critical exponents were 
found. For the solid solutions, the exponents do not depend on the concentra- 
tion of non-magnetic component and belong to the new universality class [5] 
as predicted by the Harris criterion. We do not know analogous experiments 
where an influence of disorder on criticality of easy-plane magnets was exam- 
ined. However its irrelevance was experimentally proven [7] for the superfluid 
phase transition in He 4 which belongs to the same 0(2) universality class as 
the ferromagnetic phase transition in easy-plane magnets. 

As far as the disorder should be irrelevant for the asymptotic critical behaviour 
of the Heisenberg magnets, the diluted d = 3 Heisenberg magnets should 
belong to the same 0(3) universality class as the pure ones. Theoretically 
predicted values of the isothermal magnetic susceptibility, correlation length, 
heat capacity, pair correlation function, and the order parameter asymptotic 
critical exponents in this universality class read [2]: 

7 = 1.3895 ± 0.0050, v = 0.7073 ± 0.0035, a = -0.122 ± 0.009, (2) 
rj = 0.0355 ± 0.0025, (3 = 0.3662 ± 0.0025. 

The experimental picture is more controversial. The bulk of experiments on 
critical behaviour of disordered Heisenberg-like magnets performed up to mid- 
dle 80-ies is discussed in the comprehensive reviews [8,9]. More recent exper- 
iments may be found in [10,11,12,13,14,15] and references therein. We show 
typical results of measurements of the isothermal magnetic susceptibility ef- 
fective critical exponent 7 e ff (1) in Figs. 1. As it is seen from the pictures, 
the behaviour of 7 e g is non-monotonic. The exponent differs from its value 
predicted in the asymptotic limit (2) and is a subject of a wide crossover be- 
haviour. Before reaching asymptotics 7 c g- possess maximum (except of the fig. 
l.d), the value of the maximum is system dependent: it differs for different 
magnets. 

It is standard now to rely on the renormalization group (RG) method [16] to 
get a reliable quantitative description of the behaviour in the vicinity of critical 
point. Namely in this way the cited above values (2) of the critical exponents 
of d = 3 Heisenberg model were obtained. The RG approach appeared to be 
a powerful tool to describe asymptotic [5] and effective [6] critical behaviour 
of disordered Ising-like magnets as well. The purpose of the present paper is 
to describe the crossover behaviour of disordered Heisenberg-like magnets in 
frames of the field-theoretical RG technique. In particular we want to calculate 
theoretically the isothermal magnetic susceptibility effective critical exponent 
and to explain in this way the appearance of the peak in its typical experimen- 
tal dependencies. The rest of the paper is organized as follows. In the Section 
2 we formulate the model and review main theoretical results obtained for it 
so far by means of the RG technique, effective critical behaviour is analyzed 
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in the Section 3, we end by conclusions and outlook in the Section 4. 





Fig. 1. Experimentally measured isothermal magnetic susceptibility effective crit- 
ical exponent j e g for disordered Heisenberg-like magnets (r = (T — T c )/T c ). a.: 
Fe 2 oNi 56 B24 (Fahnle et al, 1983 [11]); b.: Fe32Ni 36 Cri4Pi 2 B 6 (Kaul, 1985 [9]); c: 
Fe 2 oNi 6 oPi 4 B 6 , Fe4oNi 4 oPi 4 B 6 (Kaul, 1985 [9]); d.: Fei Ni 7 oBi 9 Sii (Kaul, 1988 [12]); 
e.: Fei 6 Ni 64 Bi 9 Sii (Kaul et al, 1994 [13]); f.: Fe 86 Co 4 Zri (Babu et al, 1997 [14]); 
g.: Fe 90 Zri (Babu et al, 1997 [14]); h.: Fe 90 - x Mn x Zrio (Perumal et al, 2001 [15]). 
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2 The model and its RG analysis 



The model of a random quenched magnet we are going to consider is described 
by the following Hamiltonian: 

h = -\Y. J (I R - R1)SrSWc rCr ,. (3) 

1 R,R' 



Here, the sum spans over all sites R of <i-dimensional hypercubic lattice, 
J(|R — R'|) is a short-range (ferro)magnetic interaction between classical 

— * — * — * 

"spins" Sr and Sr/. We consider the spins Sr to be m-component vectors 
and the Hamiltonian (3) contains their scalar product. Obviously, for the par- 
ticular case of Heisenberg spins we will put later m — 3. The randomness is 
introduced into the Hamiltonian (3) by the occupation numbers cr which are 
equal 1 if the site R is occupied by a spin and if the site is empty. Con- 
sidering the case when occupied sites are distributed without any correlation 
and fixed in certain configuration one obtains so-called uncorrelated quenched 
m-vector model. 

In principle, the above information is enough to apply the RG approach for a 
study of the critical behaviour of the model (3). One should obtain an effective 
Hamiltonian corresponding to the model under consideration and then one 
analyzes its long-distance properties by analyzing appropriate RG equations 
[16]. But already on this step there are at least two different possibilities 
to proceed and both were exploited for the model (3). On one hand, to get 
the free energy of the model one can average the logarithm of configuration- 
dependent partition function over different possible configurations of disorder 
[17]. Then, making use of the replica trick [18] one arrives to the familiar 
effective Hamiltonian [19]: 




describing in the replica limit n — > critical properties of the model (3). Here, 
fi is a bare mass, u > and v < are bare couplings and <p a = Q (R) is an 
a-replica of m-component vector field. The prevailing amount of RG studies 
of the critical behaviour of quenched m-vector model was performed on the 
base of the effective Hamiltonian (4) [5]. 

However, one more effective Hamiltonian corresponding to the model (3) is 
discussed in the literature [20,21,22,23]. It is obtained exploiting the idea that 
a quenched disordered system can be described as an equilibrium system with 
additional forces of constraints [24] . In such approach both variables 5r and cr 
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are treated equivalently and one ends up with the effective Hamiltonian which 
differs from (4) and, consequently, leads to different results for the critical 
behaviour of the model (3) [20,21,22,23]. Whereas the effective Hamiltonian 
(4) was used in the wide context of general m- vector models [5], the approach 
of Refs. [20,21,22,23] was mainly used in explanations of crossover behaviour 
in Heisenberg-like systems [25]. Below, we will discuss our results, based on 
the effective Hamiltonian (4) for m = 3 and compare them with those derived 
in [20,21,22,23]. 

As it is well known, the renormalization group (RG) approach makes use of 
the scaling symmetry of the system in the asymptotic limit to extract the 
universal content and at the same time removes divergencies which occur 
for the evaluation of the bare functions in this limit [16]. A change in the 
renormalized couplings u, v of the effective Hamiltonian (3) under the RG 
transformation is described by the flow equations: 

€^«(€) = A, («(€), «(€)), i±vW=j3 v (u(e)MQ)- (5) 

Here, £ is the flow parameter related to the distance r to the critical point. 
The fixed points (u*,v*) of the system of differential equations (5) are given 
by: 

P u (u*,v*) = 0, p v (u*,v*) = 0. (6) 

A fixed point is said to be stable if the stability matrix 

Bij = d(3 Ui /duj, i,j = 1,2; Ui = {u,v}, (7) 



possess in this point eigenvalues ui,U2 with positive real parts. In the limit 
£ — > 0, u(£) and v(£) attain the stable fixed point values u*,v*. If the stable 
fixed point is reachable from the initial conditions (let us recall that for the 
effective Hamiltonian (3) they read u > 0, v < 0) it corresponds to the critical 
point of the system. The asymptotic critical exponents values are defined by 
the fixed point values of the RG 7-functions. In particular the isothermal 
magnetic susceptibility exponent 7 is expressed in terms of the RG functions 
7^ and 7^2 describing renormalization of the field <fi and of the two-point vertex 
function with a 2 insertion correspondingly [16]: 

7 - = 1-5*1-. (8) 

2 - l4> 



In Eq. (8), the functions 7^ = 7^(u, v), 7^2 = 7^2 (m, v) are calculated in the 
stable fixed point u*,v*. In the RG scheme, the effective critical exponents 
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are calculated in the region, where couplings u(£),v(£) have not reached their 
fixed point values and depend on £. In particular for the exponent 7 c g- one 
gets: 



1 - 



% 2 [u{£(t)},v{£(t)}} 
2~ 1 4u{£(t)},v{£(t)}] 



+ .... 



(9) 



In (9) the part denoted by dots is proportional to the /3-functions (5) and 
comes from the change of the amplitude part of the susceptibility. In the 
subsequent calculations we will neglect this part, taking the contribution of 
the amplitude function to the crossover to be small [28]. 

For the effective Hamiltonian (4), the fixed point structure is well established 
[5]. It is schematically shown in Figs. 2. a, 2.b. Two qualitatively different 
scenarios are observed: for m > m c the critical behaviour of the disordered 
magnet is governed by the fixed point of the pure magnet (u* > 0, v* — 0), 
whereas for m < m c the new stable fixed point (u* > 0, v * < 0) governs 
the asymptotic critical behaviour of the disordered magnet. At the marginal 
dimensionality m c which separates these two regimes, the a exponent of the 
pure magnet equals zero in agreement with the Harris criterion. 



A 



A 




Fig. 2. Fixed points structure for the effective Hamiltonian (4) at d = 3 and 
arbitrary m. a: m > m c , b: m < m c . Stable fixed points are shown by filled boxes, 
unstable ones are shown by filled circles. Only stable fixed points with coordinates 
u* > 0, v* < are reachable for the model of the quenched magnet (3). 



Best theoretical estimates of m c definitely support m c < 2: m c = 1.942 ±0.026 
[26], m c = 1.912 ± 0.004 [27]. Consequently, the fixed point structure of the 
model of diluted Heisenberg-like magnet (m = 3) is given by Fig. 2. a: the 
stable reachable fixed points of the diluted and pure Heisenberg-like magnets 
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do coincide (u* ^ 0,v* = 0), hence their asymptotic critical exponents do 
coincide as well. However the last statement does not concern the effective 
exponents. These are defined by the running values of the couplings u(£) ^ 
0,v(£) 7^ and will be calculated in the next section. 



3 The RG flows and the effective critical behaviour 



The RG functions of the model (4) are known by now in pretty high orders of 
the perturbation theory [5,29]. For the purpose of present study we will restrict 
ourselves by the first approximation where the described crossover phenomena 
manifests itself for the Heisenberg-like disordered magnets in non-trivial way. 
Within the two loop approximation in the minimal subtraction RG scheme 
[30] the RG-functions read [31]: 



n , . . m + 8 3m + 14 2 5mn + 82 2 

Pu(u, v) — — u(e u — 2v H u H v + 

1 v ; v 6 12 36 

11m + 58 . , . 
jg— uv), (10) 

_ , . . m + 2 ran + 8 5(m + 2) 2 
p v {u,v) = -v{e —u v + — u + 

3mn+ 14 2 ll(m + 2) 

. , m + 2 2 mn + 2 2 m + 2 
^{u,v) = ^-u + ^ v +^—uv, (12) 

. m + 2 ran + 2 m + 2 2 ran + 2 2 m + 2 
lp(u,v) = — g— u+ — v — u ^— v — uv. (13) 

Here, e = 4 — d and replica limit n = is to be taken. 

Starting form the expressions (10)-(13) one can either develop the e-expansion, 
or work directly at d — 3 putting in (10), (11) e — 1 and considering renor- 
malized couplings u, v as the expansion parameters [32]. However, such RG 
perturbation theory series with several couplings are known to be asymptotic 
at best [16]. One should apply appropriate resummation technique to improve 
their convergence to get reliable numerical data on their basis. We used sev- 
eral different resummation schemes for this purpose. Here we will give the 
results obtained by the method which allowed to analyze the largest region 
in the parametric u — v space. The method was proposed in Ref. [33] and 
was successfully applied to study random d = 3 Ising model [29]. Moreover, 
it was shown that the RG functions of the d = random Ising model are 
Borel-summable by this method [33]. The main idea proposed in Ref. [33] is 
to consider resummation in variables u and v separately. Taken that the RG 
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function f(u, v) is given to the order of p loops, one first rewrites it as a power 
series in v: 

f(u,v) = j2Mu)v k . (14) 

fc=0 

Then each coefficient Ak{u) is considered as power series in u and resummed 
as a function of a single variable u thus obtaining the resummed functions 
Al es (u). Next one substitutes these functions into (14) and resums the RG 
function / in single variable v. For the resummation in a single variable one 
may use any of familiar methods. Our results are obtained by making use of 
the Pade-Borel-Leroy method [34]. 

First, applying the above described resummation procedure to the /9-functions 
(10), (11) we get the pure Heisenberg fixed point coordinates u* = 0.8956, 
v* = 0. The stability matrix (7) eigenvalues are positive at this fixed point 
(uji = 0.577, U2 = 0.147) providing its stability. Then for the resummed values 
of the asymptotic critical exponents we get [35]: 

7 = 1.382, v = 0.701, a = -0.104, 77 = 0.030, (5 = 0.361. (15) 

We do not give the confidence intervals in (15), as far as they can be estimated 
only by comparison of changes introduced by different orders of perturbation 
theory. Note however that the results (15) are in a good agreement with the 
most accurate estimates of the exponents in the 0(3) universality class (2). 
This brings about that both the considered here two-loop approximation as 
well as the chosen resummation technique give an adequate description of 
asymptotic critical phenomena. 

Before passing to the effective critical exponents let us first analyze the cor- 
rections to scaling. For the pure Heisenberg magnet, taking into account the 
leading correction to scaling results in the following formula for the isothermal 
susceptibility: 

X (r) = T T-^l + T 1 r A ), (16) 

where the correction-to-scaling exponent is given by A = ujv with uj = 
d(3 u {u)/du\ u=u * and non-universal critical amplitudes r ,Fi. For the diluted 
Heisenberg magnet the corresponding formula includes two leading corrections 
Ai, A 2 (see e.g. [12]): 

x ( r) = r / r- 7 (i + rir Al + r / 2 r Aa ), (17) 

with critical amplitudes r^T^r^. The exponents Aj are expressed in terms 
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of the stability matrix (7) eigenvalues O0i in the pure Heisenberg fixed point: 
Aj = vuji. At this fixed point, the eigenvalues of the stability matrix (7) read: 

_ d(3 u (u,v) _d(3 v (u,v) 

~ du l«V0,f*=(b ^2 — 7^ \ u *^0,v*=0- {^o) 



It is straightforward to see that the value uj\ coincides with the exponent uo of 
the pure model whereas it may be shown (see e.g. [12,31]) that the exponent 
^2 = M/V where a and u are the heat capacity and correlation length critical 
exponent of the pure Heisenberg model. On the base of the numerical values 
of the exponents (15) we get: 

Ai = 0.405, A 2 = 0.104. (19) 



Again, obtained by us in the two-loop approximation numbers (19) can be 
compared with those in the six-loop approximation making use of the data 
(2) together with the value of oo of pure 3d Heisenberg model u = 0.782 ± 
0.0013 [2]. As we have noted above, in order to get the numerical values of the 
correction-to-scaling exponents of diluted Heisenberg model it is no need to 
consider the RG functions (10)-(13) in the whole region of couplings u,v: it 
is enough to know them for the case of the pure model (i.e. for u ^ 0, v = 0). 
However, to get the effective exponents it is necessary to study complete set 
of the RG functions (10)-(13) working also in the region where both couplings 
u and v differ from zero. 

To this end we use the above described resummation technique in order to 
restore the convergence of the RG expansions in couplings u, v. First we solve 
the system of differential equations (5) and get the running values of couplings 
u(£), v(£) (10)-(13). They define the flow in the parametric space u,v and in 
the limit £ — > attain the stable fixed point value (shown by the filled box 
in Fig. 3). Character of the flow depends on the initial conditions uo,vq for 
solving the system of differential equations (5). For the model (3), the coupling 
v is proportional to variance of disorder [5] thus one can use the ratio Ivo/^ol 
to define the degree of dilution. Typical flows which are obtained for different 
ratios \v /uo\ are shown in Fig. 3 by curves 1-3. We choose the starting values 
in the region with the appropriate signs of couplings u > 0, v < near the 
origin (in the vicinity of the Gaussian fixed point u* = v * = shown by the 
filled circle in the figure). The flow No 1 is obtained for v o = 0, it corresponds to 
the pure Heisenberg model. The flow No 2 results from the small ratio Ifo/^ol 
and corresponds to the weak disorder whereas the flow No 3 is obtained for 
large |t>o/«o| an d corresponds to the stronger dilution. 

Obtained running values of coupling constant presented by flows in Fig. 3 allow 
one to get the effective critical exponents. Calculating resummed expression 
for the effective exponent 7 cff (9) along the flows 1-3 we get the results shown 
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0,0 0,2 0,4 0,6 0,8 1,0 1,2 1,4 1,6 1,8 




Fig. 3. Flows in the parametric space of couplings. The filled box denotes the stable 
fixed point u* = 0.8956, v* = 0. Curve 1 corresponds to the flow from initial values 
with vq = 0, curve 2 starts with a small ratio |vo/«o| whereas flow 3 corresponds to 
larger \v /u \. 

in the Fig. 4. Again, the curve 1 corresponds to the effective critical exponent 




1 1 1 1 1 1 1 1 1 1 1 

-25 -20 -15 -10 -5 

In / 

Fig. 4. Effective critical exponent 7 e g- versus logarithm of the flow parameter. The 
curves correspond to flows from the Fig. 3 denoted by corresponding numbers. 

of the pure Heisenberg model, whereas curves 2 and 3 provide two possi- 
ble scenarios for the effective exponents of the disordered Heisenberg model. 
Curve 2 corresponds to the weak dilution region: here, the exponent increases 
with approach to the critical point, although the crossover region is larger in 
comparison with the pure magnet (compare curves 1 and 2 in Fig. 4). This 
may lead to the peculiar situation that the asymptotic value of the exponent 
is reached earlier than the asymptotic values of the coupling. The effective 
exponents for the flows originating from non-zero ratio |t>o/wo| always attain 
the value which are larger than the asymptotic one. But the absolute value 
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of this "overshooting" for small enough |i>o/uo| is too small to be observed 
experimentally. An experimental observation of such type of 7^ behaviour 
of the disordered Heisenberg-like magnet is provided e.g. by Fig. Id. Differ- 
ent behaviour of 7 eff is demonstrated by the curve 3 in Fig. 4. Here, before 
reaching the asymptotic region the exponent possess a distinct peak. Such 
behaviour is in agreement with observed experimental data presented by Figs, 
la-lc, le-lh. The value of maximum depends on the initial values for the RG 
flows. Larger ratio |t>o/ M o| (i- e - stronger disorder) leads to the larger maximum. 
Thus, within unique approach one may explain both scenarios observed in the 
diluted Heisenberg-like magnets effective critical exponent 7 c g- behaviour. 

As we have noticed in the section 2, the crossover behaviour of random 
Heisenberg-like magnets was analyzed by means of an alternative approach 
in [20,21,22,23]. There, the quenched disordered magnet was described as an 
equilibrium one with additional forces of constraints [24]. This resulted in an 
effective Hamiltonian which differs from (4). The fixed point structure of this 
Hamiltonian differs from those given in Fig. 2 and, for different concentrations, 
leads to different crossover regimes. In particular, it predicts that there exists 
a limiting value of concentration where the critical behaviour is governed by 
Fisher-renormalized tricritical exponents [23] which coincide with those of a 
d = 3 spherical model: 7 = 2, v — 1, a — — 1, rj = 0, (3 = 1/2. There exist two 
more fixed points which may be stable in the weak dilution regime. Their sta- 
bility differs in different orders of the perturbation theory (compare [20] and 
[22]) but the numerical values of the critical exponents do not differ essentially 
at these fixed points. The maximal possible value of the effective critical expo- 
nent 7 e ff has been estimated as 7 e ff — 2.6 [22]. However, the distinct feature of 
the behaviour of 7 e fr (t) obtained in [21] is its monotonic dependence. Hence, 
the experimentally observed peaks (see Fig. 1) can not be explained within 
such approach. 



4 Conclusions 



In the present paper we used the field-theoretical RG technique to study the 
effective critical behaviour of diluted Heisenberg-like magnets. The question 
of particular interest was to explain the peak in the exponent 7 e ff as function 
of distance from T c observed in some experiments. Our two-loop calculations 
refined by the resummation of the perturbation theory series resulted in typical 
behaviour of diluted Heisenberg-like magnets 7 e g exponent represented by 
curves 2 and 3 in Fig. 4. The exponent can either reach it asymptotic value 
without demonstrating distinct maximum or it can first reach the peak and 
then cross-over to the asymptotic value from above. The strength of disorder 
is a physical reason which discriminates between these two regimes. 
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Our calculations are quite general and do not specify any particular object. In 
order to fit our curves to certain experiment one should include into consid- 
eration non-universal parameters to specify the magnetic system. The same 
concerns the flow parameter £ which as we have already noted is related to 
the distance to the critical point r. In principle such calculations may be 
done. However we want to emphasize that our analysis shows the reason of 
the peak in 7eff(/r) dependence for different disordered magnets which belong 
to the 0(3) universality class and this reason may be explained within the 
traditional RG approach. This concerns not only the magnetic susceptibility 
effective critical exponent. One more example is given by the order parameter 
effective exponents /3 c g which has minimum when r goes to zero (see e.g. [15]). 
Interpretation of this effect will be the goal of a separate study 

In conclusion we want to note that similar peculiarities of the effective critical 
behaviour may be observed in studies of disordered easy-plane magnets which 
belong to the 0(2) universality class. Since the heat capacity does not diverge 
in such systems, the RG fixed point scenario is given by the Fig. 2a as for the 
Heisenberg-like disordered magnets. Up to our knowledge such experiments 
have not been performed yet and we hope that our calculations may stimulate 
them. 

M. D. acknowledges the Ernst Mach research fellowship of the Osterreichisher 
Austauschdienst. This work was supported in part by Osterreichische Nation- 
albank Jubilaumsfonds through grant No 7694. 
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